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Abstract 

o ■ 

O 

■ The flow-equation method of continuous unitary transformations is used to 

ehminate the minimal quark-gluon interaction in the light-front quantized 



00 



X 



QCD Hamiltonian. The coupled differential equations in the two lowest Fock 



The original gauge field coupling is completely eliminated, even in the presence 



■ sectors correspond to the renormalization of the light-front gluon mass and the 

OO ■ generation of an effective quark-antiquark (as well as gluon-gluon) interaction. 

o 

OO 

o 

<^ • of degenerate states connected by this interaction. Further, a more singular 

o _ 

behavior for the quark and gluon effective interactions at small gluon 
momenta is obtained, due to the asymptotic behavior of the effective gluon 



,1^ • mass for small cutoff. We discuss the consequences of this asymptotic behavior 



and possible confinement implications. 

I. INTRODUCTION 

The perturbative aspects of Quantum Chromodynamics (QCD) were understood many 
years ago with the convincing documentation that QCD is asymptotically free. However, 
the calculational techniques for nonperturbative QCD are still under developement. Thus, 
quantitative analysis of low energy and momentum transfer phenomena remains difficult 
even though the qualitative features are reasonably described due to chiral symmetry and the 
phenomenological success of the constituent quark picture. Nevertheless, it is widely believed 
that pure Yang-Mills theory, without dynamical quarks, exhibits confinement represented 
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by a linear potential between static color sources. Confinement may be attributed to mass 
generation from transmutation of dimensions in QCD. Adding dynamical quarks also breaks 
chiral symmetry spontaneously An ultimate aim of nonperturbative QCD studies is to 
understand both confinement and chiral symmetry breaking and how the constituent quark 
picture arises from fundamental QCD. 

In the past few years there were several studies addressing the issue of confinement 
and mass gap generation in the framework of the Schrodinger picture [1,2]. Using path 
integral techniques, ref. [1] utilized a vacuum wave functional ansatz suggested by Kogan 
and Kovner, and integrated over all possible gauge configurations. Since calculations are 
quite formidable, this study only solves a field theory problem in 1 + 1 dimensions [2] and 
restricts 3 + 1 dimension [1] analyses to ground states. For earher investigations, see ref. [3] . 

Alternatively, nonperturbative studies have used a Hamiltonian framework with the QCD 
Hamiltonian quantized in a specific gauge. In particular the Coulomb gauge has recently 
been investigated [4,5]. In this paper we also use a Hamiltonian approach but utilize the 
light-front gauge, A'^ = [6]. There are arguements that the light-front gauge may be 
the most suitable framework to study nonperturbative QCD [6]. This conjecture is also 
supported by the success of phenomenological constituent light- front quark models [7] . To 
provide further insight concerning this issue, we have applied the flow equations to the light- 
front QCD Hamiltonian to generate dynamical gluons and quarks as well as their effective 
interactions. 

Previous light-front studies of confinement and bound states have been conducted us- 
ing the methods of similarity renormalization [8] and transverse lattice [9]. Significantly, 
light-front QCD in 3 -|- 1 rigorously contains a confining interaction in the form of the in- 
stantaneous four-fermion interaction, l/g"*" ^, which is the complete confining interaction in 
1-1-1 QCD for the light-front spatial dimension, x~ . Wilson suggested that a starting 
point for analyzing full QCD (with confinement) in light-front coordinates is the light-front 
infrared divergences. When an infrared (IR) cutoff is introduced, appropriate counterterms 
are necessary to restore the attending physics below the cutoff. Based on light-front power 
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counting, these counterterms can involve color charge densities and exhibit an unknown non- 
local behavior in the transverse direction which represents a possible source for transverse 
confinement. In the similarity renormalization approach, it has been claimed [8] that IR 
divergences from the instantaneous gluon exchange potential are not completely cancelled 
leading to a remnant potential that increases logarithmically with either increasing separa- 
tions \x~\ or \x±\. However, the issue of a nonzero gluon mass and local gauge invariance is 
not yet completely understood and the task of restoring rotational symmetry is still difficult 
to achieve. 

In this work, we also introduce an IR longitudinal cutoff and generate a light-front 
counterterm which sets a scale for a dynamical mass gap and string tension. While our 
calculations for the flow equation are only to order g"^, our results challenge the conven- 
tional notion that weak-coupling Hamiltonians derived from QCD have only Coulomb-like 
potentials, and deflnitely no confining interactions. Introducing a longitudinal IR cutoff in 
light-front dynamics makes it impossible to create particles from a bare vacuum by a trans- 
lationally invariant Hamiltonian and thus the number of constituents in a given eigenstate 
is fixed. However, light-front counterterms for the longitudinal IR cutoff dependence can 
generate a nonzero amplitude of particle creation and are therefore a possible source for 
features associated with a nontrivial vacuum structure in equal-time dynamics, including 
confinement and spontaneous symmetry breaking. Note that small light-front longitudinal 
momentum fraction x corresponds to high light-front energy. Therefore, in order to remove 
small X divergences and maintain the cutoff-independence, one should use a renormalization 
group, which is appropriate to high energies. 

In this work we address the above issue using the flow equations method [10]. More 
specifically, we adopt Susskind's idea of the 'long arm' of the vacuum [11] which addresses 
chiral symmetry breaking in the light-front formulation. In the parton model, one pictures 
a high momentum hadron as a collection of constituents each also having large momentum, 
such that under a boost, both the hadron and the partons change their momenta by the 
same factor. Therefore, one can formulate an effective field theory on the axis of the light- 
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front longitudinal momentum fraction x (or on the axis of rapidity which is logarithm of 
x). The partons, both valence and sea, have positive x and, according to Feynman and 
Bjorken, are distributed along the x-axis with the density of sea partons increasing for small 
X according to dx/x. Here, the vacuum is at x = 0. The fundamental property of light-front 
Hamiltonians is that under a rescaling of the light-front momentum, x ax, the light-front 
Hamiltonian scales H — > H/a. This may be interpreted as a dilatation symmetry along 
the X-axis. This symmetry only holds classically and is broken on the quantum level by an 
anomaly. Now, the long arm of the vacuum occurs because the coupling, ie. the interaction 
between neighboring partons, gets stronger and stronger as one approaches a; = so rapidly 
that the system is able to hold itself together despite the fact that there is an infinite number 
of steps between x — Q and finite x [11]. To illustrate the long arm effect of the vacuum, 
we introduce a cut-ofi^ for small x. A natural cutoff is provided by 5x = e'^x, the minimal 
spacing between constituents, which plays the role of UV-regulator. As long as the density 
of partons on the rapidity axis is not infinite, e'^ or 5x is finite and one obtains finite matrix 
elements. The cutofi^ 6x breaks dilatation symmetry in the x-axis and generates an energy 
scale, or mass gap governing the strength of the effective interaction between quarks, in our 
case a string tension in the quark-antiquark potential. Hence the long arm of the 'light-front 
vacuum' enables formation of a qq bound state through breaking of an internal symmetry, 
analogous to the creation of Cooper pairs in a superconducter. Because results must be 
independent of the cutoff, renormalization group equation is required, which in this work is 
provided by the Hamiltonian flow equations method [10]. 

Incorporating effects from small x into an effective light-front Hamiltonian is equivalent 
to integrating over the high light-front-energy modes in the asymptotically free domain. In 
terms of the renormalization group, regulating small x introduces a mass gap, which together 
with asymptotic freedom leads to a renormalization group invariant scale and dimensional 
transmutation along the x-axis. This transmutation in turn produces a linear effective 
interaction more singular than the Coulomb potential. 

At the technical level our thrust is to use the flow equations to renormalize the gluon 
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energy and eliminate the coupling between the gluon and quark Fock sectors. We focus 
upon zero gluon momentum and obtain a gap equation for the renormalized light-front 
gluon mass /x(A). This equation can be solved by imposing a renormalization condition, 
/^ren(A = Aq) = i^o, whcrc /iq IS the 'physical' value. The renormalized effective gluon 
mass, firenW, IS obtained by introducing a mass counterterm. As a result, the asymptotic 
behavior of Hreni^) is obtained for small cutoff A which approaches the renormalization 
point Ao from above (A > Aq). The properties of //ren(A) ensure the quark-gluon coupling 
is eliminated even in the degenerate case of zero gluon momenta. A similar idea was used 
for the spin-boson model by Kehrein, Mielke and Neu [12], who argued that the coupling 
to a bosonic bath is always eliminated by renormalizing the tunneling frequency. Also, in 
complete analogy to our problem, Lenz and Wegner [13] showed for interacting electrons 
in BCS-theory, that the elimination of electron-phonon coupling for all states is a direct 
consequence of renormalizing the phonon frequency. Finally, the flow of the gluon mass 
with A produces for small momentum q a potential of form 1/g^ which is more singular than 
the Coulomb l/q^. 

In the next section (Section II), we consider flow equations as a renormalization group 
transformation for Hamiltonians and formulate them for one and two-body sectors explicitly. 
In Section III the flow equations are applied to the light-front QCD Hamiltonian to generate 
a gluon gap equation and an effective qq (as well as gg) interaction. This section also 
addresses solving these equations. Concluding discussion follows in Section IV. 

II. FLOW EQUATIONS FOR GAUGE THEORIES 

A. Flow equations as renormalization group transformations 

The basic element of the renormalization group transformation is a unitary transforma- 
tion that renders the Hamiltonian matrix band-diagonal; i.e. matrix elements with energy 
differences \Ei — Ej\ exceeding a cutoff A are eliminated [14]. This procedure converges well 
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when there is a hierarchy of scales in the problem. The goal is to decouple the high- and 
low-energy scales of the band-diagonal effective Hamiltonian, which is renormalized order 
by order in perturbation theory. Using A^th order perturbation theory to connect high and 
low energy states of an effective renormalized Hamiltonian with ultraviolet regulating cutoff 
A (which is the size of the full Hamiltonian matrix in the energy space, eventually taken to 
infinity), requires N = 2(A — A)/A. Assuming a coupling constant g < 1, g'^^^~'^^^^ <^ 1 for 
A — > oo. One has therefore isolated the low energy scale effective Hamiltonian which can 
then be diagonalized nonperturbatively for the few lowest eigenstates. This should be simpler 
than solving the full Hamiltonian matrix exactly. The unitary transformation which connects 
Hamiltonian matrices with different widths, i?(A2,A) = C/(Ai, A2; A)if(Ai; A)C/'^(Ai, A2; A), 
is the renormalization group transformation formulated by Glazek and Wilson [14] called 
the similarity renormalization. 

In terms of an effective Hamiltonian, one removes large energy differences by squeezing 
the width of the band, A2 < Ai <^ A. This is done sequencially from high to low energy 
to construct a low-energy effective Hamiltonian. Since the light-front energy has an inverse 
relationship to the light-front momentum, Elp ~ IZ-P"*", this scaling to lower A is equivalent 
to moving from small x (high energy) to large x (low energy) for an effective light-front 
theory formulated on the x-axis. Elimination of large energy differences is governed by a 
dimensionless formfactor, called the similarity function, f{\Ei — Ej\/X), which is of order 
unity for small arguments (no elimination), and approaches zero for large arguments (com- 
plete elimination of the off-diagonal matrix elements). Note that squeezing the band width 
to zero is equivalent to exactly diagonalizing the Hamiltonian matrix. 

The flow equation approach entails an infinitesimal unitary transformation written in 
differential form 

r,(l) = lHi(l),HAl)\, (2.1) 
where I is the flow parameter related to the cutoff scale A by / = 1/A^. The choice for the 



transformation generator 77 has been suggested by Wegner [10] as the commutator of the 
diagonal, Hd, and off-diagonal (rest), H^, Hamiltonian components. As shown by Wegner, 
for / — > 00 the off-diagonal part of Hamiltonian is eliminated, producing an effective diagonal 
Hamiltonian. 

One of the purposes for renormalization is to remove the UV divergent high-energy 
contributions. Since they arise from the particle-number-changing terms of the Hamiltonian 
(at least in gauge theories), only these terms should be ehminated. Wegner suggested 
assigning the off-diagonal part to the particle- number-changing term Hr — Hn^m, and 
the diagonal part to the particle-number-conserving term Hd = Hn^n. Thus instead of 
diagonalization, one can implement flow equations to block-diagonalize the Hamiltonian 
in particle number space. Also, block-diagonalization in particle number space precludes 
convergence problems associated with exact diagonalization in energy space [10]. Most 
noteworthy is the practical aspect of block-diagonalization when applied to field theory. 
By block-diagonalizing the Hamiltonian in particle number space the high and low Fock 
sectors are uncoupled, enabling separate eigenstate problems with an effective Hamiltonian 
in each particle number sector. For most calculations, solving an effective sector-uncoupled 
Hamiltonian should be simpler than solving the original Hamiltonian [15]. Note, however, 
that the flow equations eliminate the particle number changing matrix elements not in one 
step, but rather continuously for different energy differences, sequencially from high to low 
energy (matrix elements between degenerate states, Ei — Ej, can also be eliminated by flow 
equations as discussed in the next subsection). Here, the link between the similarity and 
flow equation schemes and renormalization clearly emerges. The distinguishing feature for 
flow equations, however, is the separation of the particle-number-conserving and particle- 
number-changing terms, with only the particle-number-changing terms contributing in the 
renormalization of an effective Hamiltonian. 

In this work, we run flow equations for the two lowest Fock sectors {\qq) and \g)) in a 
gauge theory and obtain equations for the single gluon energy and qq effective interaction. 
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B. Flow equations in the two lowest Fock sectors 



We consider a canonical Hamiltonian operator for a gauge field theory, where the gauge 
field is minimally coupled to the matter field (for example, minimally coupled QED or simple 
abelian QCD). In terms of bare quark and gluon fields the eigenfunctions of this Hamiltonian 
contain infinitely many Fock states. For the two lowest Fock states (neglecting the rest), 
the Hamiltonian matrix in particle number space is 



H 



^ PHP PHQ^ 
y QHP QHQ J 



(2.2) 



where P and Q are projection operators for one- and two-body Fock states. For abelian, 
minimally coupled QCD, P projects on a one-gluon state and Q on a quark- antiquark pair, 
= \9) ^iid Q\i^) = \qq)- One-quark states are omitted because there are no dynamical 
quarks in this analysis (we prefer to disentangle complexities from chiral symmetry breaking) . 
Matrix elements of the operator PHQ describe quark-gluon minimal coupling, the PHP 
term is a gluon effective energy and QHQ describes the qq effective interaction. Flow 
equations for particle number changing matrix elements, hpg, are given by 

zM^l{E,{l)-E,{l)f , (2.3) 

where p and q span all (free single particle) energy states in the P and Q subspaces, re- 
spectively; Ep{l) and Eq{l) are diagonal matrix elements of block Hamiltonians PHP and 
QHQ. In the first equation the off-diagonal matrix elements hppi and hgg' are neglected to 
leading order in coupling 0{hpq) (this uncouples the flow equations for H^ and Hr terms). 
With the generator rjpg defined by Eq. (2.3), elimination of the coupling between the P and 
Q sectors 

hpg{l) = hpg{0)f{Zpg{l)) , (2.4) 



is governed by the similarity function f{l{Ep{l)—Eq{l)Y), which vanishes for matrix elements 
with energy differences exceeding the cutoff A, \Ep{l) — Eq{l)\ 3> = X. Also, hpq{0) 

is the initial value, and Ep{l) and Eq{l) flow too (see Eq. (2.6)). Eliminating the coupling 
hpq generates effective interactions in the particle number conserving P and Q sectors. The 
corresponding flow equation reads 

- E ivMKp'ii) - hpqiiM)) , (2.5) 

with a similar equation for hqq'. Using the generator defined in Eq. (2.3) these equations 
reduce to 

^-?^(0^|(""«M0) (2.6) 

dhqqjl) _ dhpqjl) 1 ^ ^^ 1 dJip,^.{l) \ 

dl ~ ^[ dl Eq{l) - Epilf''^'^'^ ^ ^'^^^'^ Eq,{l) - Ep{l) dl 

with the coupling hpg given by Eq. (2.4). Note, for the first equation in Eq. (2.6) only the 
diagonal matrix elements in the P space enter. In our application these equations represent 
the gap equation for an effective gluon energy and an equation for an effective ^'^-interaction, 
both investigated in the next section. 

Generally, the set of Hamiltonian flow equations is not finite and the equations are cou- 
pled by kernels, which also flow with that are only known after solving the flow equations. 
Obviously, practical computations require truncating to the few lowest sectors, assuming 
that higher sectors are not very important at low energies. In the light-front framework, 
this approximation is valid because pair creation from the light-front vacuum is forbidden 
and generally the higher Fock components (with large particle number) carry large light- 
front energies. Neglecting the high Fock components reduces the problem to an effective 
low-energy theory. In this way, we close the set of equations and reduce the number of 
unknown couplings, leaving only the canonical operator couplings. Note that the truncation 
in number of particles participating in intermediate states is not equivalent to perturba- 
tion theory in coupling constant, but rather is closer to the Tamm-Dancoff approach. The 
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flow equations are stiU coupled, and should be solved self consistently including the flow- 
dependences of the couphngs. 

In this work on minimally coupled gauge theory, the set of flow equations was truncated 
at the two-body sector, including at most three-particles in intermediate states. Hence 
the renormalization of the QCD strong coupling constant is not considered, but the single 
particle energy, Eq. (2.6) for Ep{l), is renormalized. In particular, as shown in the next 
section, the renormalization of a gluon energy is important for constructing an effective 
^^-interaction over all energy ranges. 

Flow equations not only eliminate the quark-gluon coupling, that renders the Hamilto- 
nian block-diagonal, they also permit renormalization. Renormalizing the gluon effective 
energy yielding Ep{l) leads to two important consequences. First, the quark-gluon coupling 
is eliminated for all energies, even for degenerate states. Second, eliminating the quark-gluon 
coupling with Ep{l) for degenerate states, the effective qq interaction is obtained of the form 
which is more singular than the perturbative result, 1/q^. This potential provides 
confinement and bound state of quarks. 

We conclude this section with a synopsis of previous flow equation applications in other 
flelds. The situation with degenerate states was flrst treated by flow equations in solid 
state physics. Lenz and Wegner showed for systems of interacting electrons [13], that the 
original electron-phonon coupling can be completely eliminated, even when the states con- 
nected by this interaction are degenerate. Lenz and Wegner found [13] that the induced 
electron-electron interaction differs from Frohlich's, whose unitary transformation is based 
on second order perturbation theory. Moreover, this interaction is attractive for any mo- 
mentum, binding electrons in Cooper pairs. Also, Kehrein, Mielke and Neu [12] have shown 
for the spin-boson problem, that flow equations allow a complete elimination of the coupling 
to the bosonic bath even for real processes. Finally, Kehrein and Mielke obtained similar 
modiflcations due to an Z-dependent generator by eliminating the hybridization term in the 
Anderson model [16]. The authors showed that their approach generates a spin-spin inter- 
action which differs from the one obtained by the famous Schrieffer- Wolff transformation. 
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Their induced interaction has the right high-energy cutoff, as compared to Schrieffer-Wolff 's 
result. Thus, within flow-equation approach it is possible to obtain new results which can 
not be obtained by perturbation theory 

In the next section, we use the above formulation to solve the flow equations for the 
effective gluon energy and the quark- antiquark interaction (see Eq. (2.6)) using the light- 
front quantized QCD Hamiltonian. 



III. FLOW EQUATIONS IN LIGHT-FRONT QCD 

A. Gluon gap equation 

In the light-front formulation, the flow equation for a single particle energy = + 
//^(A))/p"'" is actually written for the mass /x^(A) since the term p]_/p^ is independent of 
flow. The set of the coupled hght-front equations for the cutoff dependent quark and gluon 
masses was flrst derived by Glazek [17]. We uncouple this set of equations by assuming that 
the quark mass does not flow with the cutoff, m(A) — m, where m is the current quark 
mass. The light-front gluon gap equation is 



d\ ^ ^ Jo x{l-x)Jo 167r3^*' 'gi(A) dX 

"Jo x{l-x)Jo 16n^^<'^ ^Ql{X) dX 



where 



Here, the group factors for SU{Nc) are TfSab = Tr(r"r'') = ^Sab and the number of flavors 



2" 

colors Nc — 3 (the subscripts in Ca and Nc should not be confused with the group indices 



is Nf = 6. The adjoint Casimir is defined by CaSab = f"-'^'^f^^'^ = N^Sab with the number of 
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a and c). In the integral kernel, the gluon couples to quark- antiquark pairs and gluon pairs 
with quark-gluon coupling, gq{X), and three-gluon coupling, gg{\), respectively. For non-zero 
A these couplings are different functions of momenta. The light-front momentum flowing in 
the loops has components {x, k±) . In our derivation, the connection between flow parameter 
and the cutoff, I — is used. 

In Eq. (3.1), the effective gluon mass is defined at transverse gluon momentum = 0, as 
in light-front perturbation theory for the gluon mass correction [18]. Following other gluon 
gap equation studies, we assume the effective gluon mass vanishes for large gluon momenta. 
Therefore, even though the effective gluon mass depends on momentum, only its A cutoff 
dependence at zero momentum is given by Eq. (3.1), i.e. lJ'{X,p{Pz,P±) — 0). However, the 
A dependence of /j, is the only relevant renormalization issue. 

Generally, it is difficult to solve Eq. (3.1) because the running couplings, gq{X) and 
gg{X), are not known. Therefore, the couphng cutoff dependence is neglected. Also, the 
initial condition for Eq. (3.1) is not known. Accordingly, the following two renormalization 
conditions are imposed to determine the running gluon mass n{X). First, the effective 
Hamiltonian at scale A has eigenstates with eigenvalues p~ — {p']_ + /J'D/p'^ [17], satisfying 

^^(p \p) = — ^ — {p \p) 

- /'ciAX(V.(A')V(A') - V.(A')^.p(A')) , (3.3) 
q 

where /^o denotes the 'physical' gluon mass, and \p) denotes a single effective gluon state (P 
subspace, Eq. (2.2)) with momentum {p'^,p±) and {p'\p) — 167r^p+5^^)(p'— p). The generator 
rjpq, given by Eq. (2.3), eliminates the quark-gluon (three-gluon) coupling hpg. Second, the 
effective gluon mass, renormalized in second order perturbation theory, equals the 'physical' 
mass 

//L(A = Ao)=//^, (3.4) 

at the renormalization point Aq. 
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Prom the definition of fiow equations, a resulting effective Hamiltonian is given in the 
hmit of the flow parameter / ^ oo, or A = 0. In the parton picture described in the 
introduction, the boost may be regarded as a renormahzation group operation and the 
renormahzation group flxed point can be identifled as the inflnite momentum hmit. The 
corresponding renormahzation group flxed point in the flow equations is then A = 0. There- 
fore, in Eq. (3.4), the renormahzation point is set to zero, Aq = 0, at the end of calculations. 
Also, the 'physical' gluon mass, which is used as an arbitrary mass parameter, can be taken 
to zero to restore gauge invariance (see below). 

Even with these simpliflcations, solving Eq. (3.1) is still quite involved. The solution can 
only be found numerically. In order to proceed analytically, a mass parameter /iq ('physical' 
mass), that corresponds to the choice of the renormahzation point at = ^l, is substituted 
for /^(A) in the integral kernel on the right-hand-side of Eq. (3.1). The same procedure is 
used to calculate the perturbative mass correction in [19]. Therefore the energy differences 
in Eq. (3.2) are given by 

Oh4^-^o, ^^ = Tr^-^o- (3.5) 

x{l — X) x{l — X) 

Using Eq. (3.5), Eq. (3.1) is solved iteratively and A*^(A) — /Xg is the flrst iteration. 

Another problem in Eq. (3.1) is that the loop integrals have an UV divergence in the 
transverse directions k± and an IR divergence in the longitudinal direction x. The flow 
equations naturally regulate these divergences via the similarity function in each vertex (for 
example, the gluon loop in Eq. (3.1) is regulated by Qf(A) < A^). This type of regularization 
is known as 'Jacobi cutoffs', because of the Jacobi momenta of a constituent (it is also 
called 'global' regularization in [19]). The advantage of this regularization in the light- front 
approach is that it preserves both transverse and longitudinal boost invariance [18]. For a 
nonzero mass such as a current quark mass in the quark loop, this regularization restricts 
k± and the light-front x integrations. Thus, the similarity function f{Ql; A), with Ql given 
by Eq. (3.5), ensures that < k± < {X^ + iJ'l)x{l — x) — w? and m^/(A^ + /^o) ^ 2; < 
1 — m^/(A^ -|- //q). However, for zero gluon mass in the gluon loop, k±_ is restricted from 
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above and x runs the entire range < x < 1, contributing divergent terms at a; = and 



included, the gluon loop is still divergent in x [19]. For a massive case the cutoff A has a 
lower bound from a mass in the theory, A > m 7^ 0, that limits the x-integration. The reason 
why flow equations do not regulate the light-front divergences in a massless case is because 
the cutoff A can evolve to A = with the similarity function f[z —>■ 00) = 1, and x spans 
the entire to 1 range, where the loop integral with massless intermediate states diverges at 
X — and x — 1. To regulate these divergences in a massless case, Zhang and Harindranath 
suggested restricting the k± integration also from below by some scale u. This is equivalent 
in our approach to integrating the flow equation, Eq. (3.1), for finite limits, from A to w [19], 
which restricts intervals in the x-axis to above zero and below one. It was also shown in 
[18] and [19] that for the scale u, even if /xq = is zero, the gluon mass correction does not 
vanish. Effectively, introducing the scale u mimics the situation of a non-zero mass, m — u, 
in intermediate states. 

In terms of effective light-front theory formulated on the x-axis, the Hamiltonian below 
the first light- front cutoff, i?o<x<£, as well as the Hamiltonian above the second cutoff x — 1, 
i7(i_e)<a;<i, describe the strongly correlated high energy states. Thus, they can be replaced 
by the Hamiltonian vacuum expectation value, since strongly coupled configurations are 
frozen. This v.e.v. replacement in gluodynamics is consistent with a composite field, 0, 
creating 0+ glueballs having finite v.e.v. [2]. The Hamiltonian in the intermediate region, 
He<x<{i-e), is then treated by flow equations. 

Integrating the gluon flow equation, Eq. (3.1), for finite limits [u; A], yields 




Even when the instantaneous diagrams are 
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y kj_- x{l - x)hq + fjLo k'i- x{l - x)iil + 111) 

where the renormahzation point is = fi^, and the energy differences are defined in 
Eq. (3.5). Though the quark loop does not diverge with x, the gluon flow equation is 
integrated the same way for the gluon and quark loops. Similarity functions restrict the 
transverse momenta to k±min ^ k± < k±max, with 

k±min = {U^ + fJ-Dxil -X)- 111, (3.7) 

and for x 

^^<x<l-J^. (3.8) 

Analagous expressions hold for the quark loop. 

Integration over momenta {x, k±) in Eq. (3.6) yields 

/.^(A)-/.>)= -^^(l/x^ln^ + mMn^ + l(A^-^.^)) (3.9) 

Assuming the renormahation conditions of Eqs. (3.3) and (3.4), the effective gluon mass at 
scale A is 

li\X)=lil + 6hIt ( A) + (5/^2 ( a, Ao) . (3.10) 

Here, Eq. (3.3) fixes the integration constant (when integrating the fiow equation (3.1)) to 
/Xq. In Eq. (3.10) the perturbative term reproduces the result from light-front perturbation 
theory [19], i.e. 
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Renormalizing the effective Hamiltonian to second oder 0{g^), the perturbative mass cor- 
rection is absorbed by the mass counterterm, m^rpi^Kuv) = — 5/ip2.(A[/y) with Aj/y oo, 
and the renormahzed effective gluon mass is 

/^L(A)=/^'(A)+m2,^(A), (3.12) 

for A = Kijy oo. Though perturbative renormahzation is apphed at large cutoff scales, 
Ajyy, we assume that the leading cutoff dependence in second order is absorbed by the mass 
counterterm for all A, and the renormahzed gluon mass is given by Eq. (3.12) for any A. 
Explicitly from Eq. (3.9), the renormahzed gluon mass reads 

f4enW = /^O + <^/^^(A> Ao) = /^o + ^(/^O, u) ^ T2 (3.13) 

where the renormahzation condition, Eq. (3.4), at scale Aq is satisfied. Recall, this solution, 
Eq. (3.13), describes an effective gluon mass at zero gluon momentum. Since the effective 
Hamiltonian generated by flow equations is defined at cutoff scale A — > 0, the resulting gluon 
mass equals the 'physical' mass, lJ>fen{X = Aq = 0) = //q- In particular, when the 'physical' 
mass is set to zero, Hq — 0, the effective QCD Hamiltonian has zero mass gauge fields, and 
our unitary transformation becomes gauge invariant. Then the Kinoshita-Lee-Nauenberg's 
theorem [23] is applicable, according to which IR singularities associated with soft collinear 
gluons are cancelled exactly only for color singlet, gauge invariant matrix elements. This was 
demonstrated in detail perturbatively for e~^e~ annihilation amplitude [24]. In our situation 
this zero gluon mass, gauge invariant limit also produces IR divergent quark and gluon self- 
energies. This in turn generates color confinement because the infinite self-energies suppress 
the quark and gluon propagation amplitudes. However, and, consistent with the above 
theorem, in the Bethe-Salpeter equation (and the simpler Tamm-Dancoff equation) the IR 
divergent contributions from the kinetic self-energy and potential interaction kernel parts 
cancel exactly for color singlet states. Our generated two-body qq and gg potentials satisfy 
these conditions providing finite masses for singlet states of mesons and glueballs. 
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Prom Eq. (3.13), one has in the limit //q — > 



a — lim (t(//o, u) 




(3.14) 



27r2 



where, as shown below, a plays the role of the string tension between quark and antiquark. 
In Eq. (3.13), v? determines the rate at which the effective gluon mass asymptotically 
approaches the 'physical' value, ~ tt^ln(A/Ao). 

Though this paper addresses confinement, we are not yet able to suggest a physical 
picture for how confinement emerges in the light-front frame. However, we note that QCD 
has a non-abelian gauge group and strong interactions, which distinguishes it from QED. 
Both canonical light-front QCD and QED Hamiltonians have the instantaneous interaction, 
1/x^, which has a singular behavior for small x. For QED a small x regularization scheme 
can be found where divergent contributions from this instantaneous term are canceled in 
the matrix elements [19]. This is not true for QCD [19]. Due to the non-abelian triple-gluon 
vertex it is necessary to introduce a regulating scale u for the divergences at x ~ and 
a; ~ 1 in the gluon self-energy (the quark loop is regulated in the same fashion in order to 
combine quark and gluon loops). In fact, the string tension in Eq. (3.14) is proportional 
to the Casimir operator in the adjoint representation refiecting this non-abelian character. 
There is no such scale in the QED self-energy operators, which are all regulated by a current 
electron mass m. One can therefore argue that the scale u carries information about strongly 
interacting fields and is present only in QCD. 

In asymptotic free theories, such as QCD, the regulating scale can be related to the gauge 
invariant scale using the Callan-Symanzik equation. In our case, the scale u can be expressed 
in terms of h.QCD by solving the third order fiow equations for the strong coupling constant, 
as{\)- Such calculations have been recently performed for an asymptotic free toy model 
[20] and for the three-gluon vertex in QCD [21]. However, here we forgo this procedure and 
simply assume u is given by the hadron scale, u ~ ^hadrom that also establishes the string 
tension scale. 

The asymptotic behavior of the effective gluon mass near the renormalization point is 
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important for selfconsistently solving the flow equations for effective interactions, Eq. (2.6), 
at vanishing gluon momenta. In the next subsection, this dependence (Eq. (3.13)) is used 
to find an effective quark potential generated by the flow equations. 

B. Effective quark-antiquark interaction 

Eliminating the quark-gluon couphng generates an effective interaction in the quark- 
antiquark sector (the second equation in Eq. (2.6)). The energy transfers (i.e. energy 
denominators) along the quark pi — > p2 and antiquark ~^ Pi hues are given by 
Di = pi — P2 — {pi — P2)~ and D2 = P2 — p'l — {p'2 ~ P'l)^ 1 respectively. We denote 
the exchange momentum as q = {q'^,q±) = pi — P2- These energy transfers can be re- 
lated to the corresponding square of the four-momenta as Q\ — (pi — p-if — —q'^Di and 
Ql = iP2 ~ P'lf = —(1^D2 [22]. In the hght-front formulation, the scalar product is deflned 
p ■ k = {p'^k~ + p~k~^)/2 — p± ■ k± and the four-momentum transfers for quark momentum 
{x, k±) and antiquark momentum (1 — x, —k±) are given by 

Ql{X)^Ql + l^lenW 

Ql{X) = Ql + l^lenW , (3.15) 

where the asymptotic form for the rcnormahzed effective gluon mass, /i^g„(A) = A*o + 
a{iio, u) InlX^ /Xq) (Eq. (3.13)), enters. The four-momentum transfers are now dependent 
on the cutoff A and for zero gluon mass are given [22] by 

2 {x'k± — xk'j^"^ + rn^lx — x')"^ 

Qi — ; 

xx' 

Ql = Qi\x^{i-xy, x'^{i-x') ■ (3.16) 

As shown below, in the instant formulation, energy transfers with zero gluon mass, 
Eq. (3.16), describe a three-momentum exchange, q^ with q{qz,(l±), while energy trans- 
fers with non zero cutoffs, Eq. (3.15), correspond to an effective energy (with effective 
mass) exchange, q^ + /x^g„(A). Near the renormalization point the energy transfers given 
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by Eq. (3.15) and Eq. (3.16) coincide asymptotically, i.e. for A — > Aq ~ and /lo ~ 0, 
Qi{X) — > Qf, i = 1,2. Using Eq. (3.15), the gluon four-momentum can be written as 

where the light-front unit vector ry^ is rj'^ = (2~, 0+, 0^), with r) ■ k — k'^. In subsequent 
calculations, we use the average and difference of four-momenta transfers 

g^(A) = iQliX) + Ql{X))/2 = Q' + /.^(A) 

SQ'iX) = (Q?(A) - Ql{\))/2 = 5Q' , (3.18) 

where = (Q? + Ql)/2 and SQ^ = (Q? - Ql)/2. 

Solving the flow equations selfconsistently, the energy transfers for non-zero A define 
effective interactions between quarks as specified by Eq. (2.6). Thus, the similarity formfac- 
tors in each vertex f{Ql{X); A) and f{Ql{X); A) provide an effective interaction (see below). 
However, for large A only high energies are eliminated and the momenta transfer cutoff 
dependence, Q^{X), is minimal. Reducing A, energy transfers are eliminated continuously 
from high to low energy. Only for small A, and corresponding momentum transfers, is the 
asymptotic behavior of Qf{X) essential. This is important in solving for the effective in- 
teraction, since (see subsection 2.2) this asymptotic behavior (the cutoff dependent gluon 
effective energy at small A) eliminates the quark-gluon coupling even for vanishing gluon 
momenta, i.e. when Ql = Ql = {x = x' and k± = fc^). If = and /xq — 0, the 
similarity formfactor / = exp(— (5^(A)/A^), with Q^{X) given by Eq. (3.15), decays at small 
A — > (with A > Ao) as 

/-^y . (3.19) 

Note this is a power suppression, not exponential, with a specified by Eq. (3.14). Only for 

A = Aq does / = 1, however, for other values, A 7^ Aq, / decays and the integral over A in 
Eq. (2.6) is finite in the effective interaction. 
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In the light-front formulation, an effective quark-antiquark interaction is given by a ma- 
trix element of an effective light-front QCD Hamiltonian in the gg-sector (Q-space) and in- 
cludes dynamical interactions generated by flow equations (the second equation in Eq. (2.6)) 
as well as instantaneous terms present in the the light-front gauge. These terms are embod- 
died in the effective g^-interaction 

V,, = -47ra,C/(7 V) , Um ^ B,,, , (3.20) 

where C/ = T"T" = (A*"^^ — l)/2Nc is Casimir operator in the fundamental representation, 
and the current- current term in the exchange channel is given by (quark helicity notation is 
suppressed) 

{u{-k^, (1 - x))^^u{k^, x)) [v{k'^, x')r'v{-k'^, (1 - x'))) 

(tV ) = / , w ^ ■ (3.21) 

lxx'[l — x)[l — x') 



The interaction kernel, B^^i, is analogous to the effective electron-positron interaction in 
the light-front QED [22], except for the cutoff dependence in the energy transfers 

B^^' = g^^' (/i + h) + ^mV^ (^1 - . (3.22) 

Here gf^^i/ is the light-front metric tensor and is given previously (see Eq. (3.17)). The cut- 
off dependent energy (four-momentum) transfer, Eq. (3.15), along the quark and antiquark 
lines is included in the integral 

^-/:-^^^/(«^(^)^^)- (-3) 

I2 is obtained by interchanging indices 1 and 2 in Eq. (3.23). As discussed previously, the 
sensitivity to the renormalization scheme and renormalization point is eliminated by setting 
both /iq and Aq to zero after integration. 

For the following three explicit forms for the similarity function, the leading behavior of 
the integral factor, Eq. (3.23), is given by 

(1) Exponetial 
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1 ( , (^{f^o,u) 



(2) Gaussian 

/ = exp(-Q(A)/A), ^. = g4(,j + Q4(,„) (l+Qf(A^ 

(3) Sharp 



f^e^X'-Q^m h = (l + ^) . (3.24) 

where from Eq. (3.15) the four-momenta are given by Ql{\o) — Q^+//o with i — Several 
properties of the similarity function were used to approximate the integral in Eq. (3.23). The 
similarity function f{z) decays for arguments z > 1, thus the integral saturates for values 
< < ((5?(Ao),Qi(Ao)) at small energy transfer. Also, //o ~ Ao ~ 0. Note Eq. (3.24) 
has the form consistent with the exchange of a gluon with an effective mass parameter 
/io between quark and antiquark at distances r ~ l/z^o- The resulting gg-interaction is 
defined in the limit of vanishing mass parameter //q — > and at the zero renormalization 
cutoff point Aq — > 0, Eq. (3.20). In this limit, the average four-momentum transfer reduces to 
hm(^o,Ao)-»o <3^(Ao) = Q^, Eq. (3.18), and the string tension is given by hm(^g,Ao)^o (^(//o, u) = 
(7, Eq. (3.14). Using Eqs. (3.22) and (3.24), the resulting interaction kernel for the three 
similarity function choices is 

(1) Exponential 
hm 5^^/ = fM/^' T;^ + Til + 



(McAoHo ""^^ QV \Q'' J Q^Q'-5Q^ 

(2) Gaussian 



(3) Sharp 

where = (Q^)^ and ^Q*^ = (5Q^)^. The four-momenta in Eq. (3.25) can be represented 
in the 'mixed' light-front and instant representations as 
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Q2 ^ ^2 --{2x- l){2x' - l)(Mi - M^f 
^ -\{x - x'){Ml - Ml) , 



(3.26) 



where, together with the hght- front momentum (x, /cj,), the instant momentum k — {kz,k±) 



enters with the connection 



X = 




\Jk'^-\- m? 



(3.27) 




+ kl- The gluon three-momentum transfer is given hj q = k — k' — (qz, q±) 



and the total energies of the inital and final gg-states are given by 



(3.28) 



Our results, Eq. (3.25), are not rotationally invariant due to truncation of Fock-space. 
However, on the energy shell (Mf = M| and SQ"^ = 0), the second term in Eq. (3.25) is 
identically zero and the effective qq interaction is given by the first term 



It is significant to note that all three choices yield the same effective interaction. Here 
(7/ = aasCjjl with a given by Eq. (3.14). Further, and interestingly, we also obtain a 
Cornell type interaction with Coulomb and linear confining potentials 



If the external particles are all on the energy shell, our light-front formulation yields a ro- 
tationally invariant potential as expected. The confining term in Eq. (3.30) (or Eq. (3.29)), 
arises from elimination of the quark-gluon coupling at small gluon momenta, which is gov- 
erned by the asymptotic behavior of the effective gluon mass at small cutoff. It is also 
significant to note that the same quark effective interaction, Eqs. (3.29,3.30), emerges in the 
limit of static (infinite heavy m — > oo) quarks. 




(3.29) 




(3.30) 
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The Nc behavior of the confining term in Eq. (3.30) with the string tension given by 

_ g'C^Cf _ g\Nl - 1) 

can be compared with available 2 + 1 lattice data [25] for the expectation value of the 
Wilson loop in the fundamental representation, (Wf{C)) = exp(— cr2+iS'c), where Sc is 
the area of the loop C . Though the data are for 2 + 1 the same Nc behavior is ex- 
pected for 3 + 1 and higher dimensions [25] . Monte Carlo calculations of the string ten- 
sion in 2 + 1 give the values ^/a^/g'^ = 0.335, 0.553, 0.758, 0.966 for the gauge groups 
SU{2), SU{3), SU{4:) and SU{5) respectively [25] (note that the coupling constant in 
2-1-1 has dimension, g"^ ~energy). The corresponding values calculated from Eq. (3.31) are 
2Ti^/{ug'^) = 0.345, 0.564, 0.772, 0.977. We see that there is excellent agreement (up to 
~ 3 %) for this Nc behavior between Eq. (3.31) and the Monte Carlo results. It is further 
interesting to notice that our analytic expression for the string tension, Eq. (3.31), has the 
appropriate Nc dependence as expected from large Nc calculations. 

Similarly the quark-quark potential can be obtained and is related to the quark-antiquark 
potential by 



= (3-32) 



since only the commutator [b^ba, b^ba^] contributes to Vgg, while [b^ba, d^da^] and [d^da, b^ba^] 
contribute to Vqq. Assuming the relation Eq. (3.32), Basdevant and Boukraa [26] showed 
that the ground state of baryons can be calculated with good accuracy. 

Note, because we use perturbation theory, we do not claim to have completed the deriva- 
tion of a confining potential in QCD. However, our 1/g ^ quark-antiquark potential might 
be a precursor to quark confinement. Perhaps equally as important, our results confiict with 
the notion that a weak-coupling expansion will never produce a confining potential. 

Finally, there are corrections O {SQ"^ / Q^) to the leading effective interaction, which de- 
pend on the direction q approaches zero (as was investigated in [22]). An important limiting 
case is the coUinear limit 
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q+ ^0. (3.33) 

This is special for light-front calculations and may cause divergences. Prom Eq. (3.26) in 
this limit, 5Q^ ~ q'^, and for sufficiently smooth similarity functions f{z), like exponential or 
Gaussian, the effective interaction Eq. (3.25) does not contain a coUinear singularity, because 
(VhVh'/q^^) ■ is finite. Thus the interaction only becomes singular if q approaches zero as 
1/q^ ('Coulomb' singularity) or l/q"^ ('confining' singularity), namely the leading behavior 
given by Eq. (3.29). Note that for mass spectrums both singularities are controllable. The 
'Coulomb' singularity is integrable, because the integral / d^q/q^ is finite for small q and 
the 'confining' singularity can be regulated in the IR. However, this is not true for the 
sharp cutoff where the rj^rini term diverges as because it appears in the combination 

{VijXIh' / Q^"^) ■ \^Q'^\- We do not attach physical significance to this singularity, which is 
a consequence of the arbitrary cutoff choice leading to a singular unitary transformation 
generator at small momenta Eq. (2.3). In summary, the coUinear singularity is completely 
absent for a smooth cutoff and only the rotationally invariant part of the effective interaction 
Eq. (3.29) remains in the limit g — > 0, where the dominant contribution to the spectrum 
is obtained. Corrections from the rj^rj^i term are finite and shift the entire spectrum by a 
constant since they are diagonal in spin space. 

As numerically documented in the positronium QED effective interaction application 
[22] , rotational symmetry holds with high accuracy even if the on-energy-shell condition for 
the external particles is removed. This holds for smooth cutoff functions and even for a 
sharp cutoff if the coUinear singular part is subtracted. Based on our analysis here, we may 
anticipate similar results for the QCD effective interaction. 

C. Effective gluon-gluon interaction 

In full analogy with section 3.2, eliminating the three-gluon coupling generates an effec- 
tive interaction between two gluons. Again, in leading order, the asymptotic behavior of the 
effective mass of the exchanged gluon is included. Using the same notation as section 3.2, 
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an effective gluon-gluon interaction is again generated by tfie fiow equations. It contains 
botfi dynamical and instantaneous gluon excfiange and is given by 

y,, = -47ra,C„(r^r^V lim B^^,, (3.34) 

(mo,Ao)^0 

wfiere now CaSab — f°^f^'^ — NcSab is tfie Casimir operator in tlie adjoint representation. 
Here, tlie current-current term in tlie excfiange channef is given by 

r^r^ = > , == Pi epU^2^^' U^i V , 3.35 

where 

rM-P(_^,p^, -P2)r^''^'^'(g,p'i, -P^) = ((Pi +v2Tr + {-'i'Vxrr + (-2P2) V) (3.36) 

X {{A^^^r'g''''' ^{-y^Yg^'""' ^{-y^YQ""''') , 

and tfie figfit-front momenta are pi — {x,k±), p[ — {1 — x, —k±) and p2 — {x',k'j_), P2 — 
(l — x',—k'j_). The gfuon pofarizations are omitted for simplicity. The interaction kernel B^^/ 
is given by Eq. (3.22), which is again defined for different cut-off functions in Eq. (3.25). The 
polarization vector property, q ■ e — r) ■ e — 0, and representation for the gluon momentum, 
Eq. (3.17), have been used to simplify the above expressions. 

Both quark- antiquark and gluon-gluon interactions have the same kernel Bi^^i, and differ 
only in prefactors. Following arguments detailed in section 3.2, the leading gluon-gluon 
interaction is also independent of cut-off function form and is given by 

Vg, = (m^) (-C„^ + (J„ • r) , (3.37) 

which is, again, a Coulomb plus linear potential between two color sources in the adjoint 
representation with aa = aasCa/2. The only difference between the qq, Eq. (3.30), and gg, 
Eq. (3.37), potentials is given by the ratio of the fundamental to adjoint Casimir operators 

V,,/V,, = Cf/Ca = (AT^ - 1)/27V2 = 4/9 , (3.38) 

for SU{3), since both potentials have the same Coulomb plus confining behavior. 
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IV. CONCLUDING DISCUSSION 



An effective QCD Hamiltonian in tfie ligfit-front gauge fias been obtained, solving tfie 
flow equations for the two lowest Fock sectors selfconsistently. It has been shown that it 
is possible to eliminate the minimal quark-gluon interaction by using a continuous unitary 
transformation. In this elimination, the coupling functions of the Hamiltonian described 
by the flow equations are renormahzed. In the two lowest Fock sectors this change of 
the couplings corresponds to the renormalization of the one-particle energies and to the 
generation of effective interactions between gluons and quarks, in particular the quark- 
antiquark interaction. In obtaining the flow difTerential equations, intermediate states with 
more than three particles were omitted. Fock number truncation in intermediate states is 
quite different from a perturbation coupling constant treatment and more similar to the 
Tamm-DancofF approach. 

Our approach has several advantages. First, the original gauge field coupling can be 
completely eliminated, even when the states connected by this interaction are degenerate. 
The continuous transformation is designed such that the transformed Hamiltonian does not 
contain any interactions between one (anti)quark and the creation or annihilation of one 
gluon. These unwanted interactions, connecting states with energy differences less than a 
cutoff scale \Ep — Eq\ < A, are present in the similarity renormalization approach because 
single particle energies are not renormahzed. They can mix low and high Fock sectors, and 
are not amendable to a perturbative treatment. Further, ignoring these low-energy interac- 
tions may break gauge invariance, and, in the light-front formulation, rotational symmetry 
as well. Second, an effective quark-antiquark interaction is rotationally symmetric when the 
external particles are on the energy shell. At small gluon momenta q, also the coUinear 
singular terms ~ l/g"*" and ~ l/g"*"^ cancel. Most interestingly, in addition to a Coulomb 
term which can be obtained in second order, there is also a more singular, confining 
term 1/q'^. Our induced ^^-interaction also differs from the similarity renormalization result, 
where the coUinear singular part of the remnant instantaneous interaction (~ l/^"*"^) yields 
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a logarithmic potential which is not rotationally symmetric. These differences stem from 
the coupling's flow parameter dependence. 

Further, due to complete elimination of the quark-gluon coupling, the flow equation for 
an effective quark-antiquark interaction can be integrated for all cutoffs including A = 0. 
In the similarity renormalization approach one removes couplings perturbatively reaching 
a minimum scale (cut-off), below which perturbation theory breaks down. The value of 
this cutoff depends on the problem considered, and might be ambiguous. For QCD this 
cutoff introduces a scale in the theory, which breaks gauge and rotational invariance [8] . In 
our approach, the regulator of small light-front x establishes a scale in the effective theory 
corresponding to the string tension in the effective quark-antiquark potential. Besides the 
nonzero scale, the resulting renormalized gluon mass vanishes asymptotically, maintaining 
gauge invariance. Also, the effective quark-antiquark interaction is rotationally symmetric 
when the external particles are on the energy shell. A shortcomming in this work is that the 
small light- front x cutoff scale u enters as an input parameter, fitted to the string tension from 
lattice calculations. This can be improved by relating the cutoff u to the renormalization 
group invariant scale Aqcd which necessitates higher Fock sector intermediate states. This 
would also permit confirmation of our regularization scheme independent results. 

Our ultimate goal is to solve the coupled chain of flow equations in different sectors 
self consistently. As shown in this work, even an approximate solution of the gluon gap 
equation combined with the flow equation for the effective qq interaction provides new 
information beyond standard perturbation theory. The next step is to address the quark 
sector, formulating the quark gap equation and obtaining the renormalized light-front quark 
mass and attending improved effective quark interaction. 

Finally, reflecting the above discussions, we note that the light-front flow equations avail 
themselves to selfconsistent computer generated solution. It would be quite interesting to 
confront hadronic data with large scale, but feasible applications of this approach. 

In summary, the light-front formulation appears quite viable for examining the long range 
aspects of QCD. Implementing the flow equations within this framework has generated new 
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insight regarding nonperturbative phenomena including confinement. 
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